Processing math: 100%

cholesky_rank1_update Interface

public interface cholesky_rank1_update

Module Procedures

private subroutine cholesky_rank1_update_dbl(r, u, work, err)

Computes the rank 1 update to a Cholesky factored matrix $A = R^T R$ such that $A_1 = A + \vec{u} \vec{u}^T$ .

Arguments

Type	Intent	Optional	Attributes		Name
real(kind=real64),	intent(inout),		dimension(:,:)	::	r	On input, the N-by-N upper triangular matrix $R$ . On output, the updated matrix $R_1$ .
real(kind=real64),	intent(inout),		dimension(:)	::	u	On input, the N-element vector $\vec{u}$ . On output, the rotation sines used to transform $R$ to $R_1$ .
real(kind=real64),	intent(out),	optional,	target, dimension(:)	::	work	An optional argument that if supplied prevents local memory allocation. If provided, the array must have at least N elements. Additionally, this workspace array is used to contain the rotation cosines used to transform $R$ to $R_1$ .
class(errors),	intent(inout),	optional,	target	::	err	An error object to report any errors that occur.

private subroutine cholesky_rank1_update_cmplx(r, u, work, err)

Computes the rank 1 update to a Cholesky factored matrix $A = R^H R$ such that $A_1 = A + \vec{u} \vec{u}^H$ .

Arguments

Type	Intent	Optional	Attributes		Name
complex(kind=real64),	intent(inout),		dimension(:,:)	::	r	On input, the N-by-N upper triangular matrix $R$ . On output, the updated matrix $R_1$ .
complex(kind=real64),	intent(inout),		dimension(:)	::	u	On input, the N-element vector $\vec{u}$ . On output, the rotation sines used to transform $R$ to $R_1$ .
real(kind=real64),	intent(out),	optional,	target, dimension(:)	::	work	An optional argument that if supplied prevents local memory allocation. If provided, the array must have at least N elements. Additionally, this workspace array is used to contain the rotation cosines used to transform $R$ to $R_1$ .
class(errors),	intent(inout),	optional,	target	::	err	An error object to report any errors that occur.